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interplay of probability theory, group theory and quantum theory. 
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1 Introduction 



It is well known that the photon number distribution of the ordinary coherent states ^, Q 
is the Poisson distribution, one of the most fundamental probability distributions, which 
governs random events (such as radioactive decays) occurring in a time (space) interval. As 
we will show in this paper the relationship between the coherent states in quantum optics 
and the probability distributions are neither coincidental nor superficial but essential. The 
main purpose of the present paper is to give unified probabilistic interpretation of the various 
coherent states. 

For the elementary binomial distribution of the probability theory, corresponding to the 
binomial expansion (1 + x)^^ = J2o^ (^)'^"' have 5^(2) coherent states (the 'binomial 
states' (BS) 0) based on the spin M/2 representation. For the multinomial distributions 
corresponding to the multinomial expansion 

(i+xi + ---+x.)^= E — -xr-.-x-, (1.1) 

no+ni-\ hnr=M ^ ' 

we have certain types of su{r + 1) coherent states. These coherent states are known for 
some time but the probabilistic interpretation seems new. Let us call them multinomial 
states (MS). They are based on the symmetric representations corresponding to the Young 
diagram 

rrn ---m m boxes. (1.2) 

In probability theory the non-compact version of the binomial distribution is well known 
and called negative binomial distribution. In this paper the negative binomial states (NBS) 
of quantised radiation field will be introduced in a parallel way as the binomial states. It 
will be shown that they are the well known coherent states of su{l, 1) algebra 0, ^ 0, ||, the 
non-compact counterpart of the compact su{2) algebra. They belong to the discrete series of 
irreducible representations. Similarly the negative multinomial states (NMS), the coherent 
states of su{r, 1) algebra belonging to discrete symmetric representations, will be introduced 
in terms of the negative multinomial distributions. It is easy to see that in certain limits 
these coherent states reduce to the ordinary coherent states and their tensor products, since 
the (negative) binomial and (negative) multinomial distributions tend to the Poisson and 
multiple Poisson distributions. 

This paper is organised as follows: In section 2 the negative binomial states are intro- 
duced directly as a square root of the negative binomial distribution. In other words they 
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are constructed in such a way that their photon number distribution is the negative bino- 
mial distribution. Then these coherent states are shown to have the displacement operator 
forms. Namely, they are created by the action of the unitary operators in SU{1, 1) acting on 
certain highest (lowest) weight states ("vacuum"). In section 3 we relate the inhomogeneous 
representation of su{l, 1) suggested by the negative binomial states to the symmetric two 
boson realisation. The two boson formulation provides natural interpretation and more ex- 
plicit formulas than those of the formal representation theory of sm(1, 1). At the same time 
this section uncovers some Lie algebraic structures hidden in the probability distribution. 
The physical and statistical properties of the NBS as well as their dynamical generation are 
discussed in some detail in our recent publication 0. Section 4 deals with the generalisation 
to su{r,l), the negative multinomial states. One formulation of the negative multinomial 



states is closely related with the Holstein-Primakoff (H-P) [|T0| type realisation of su{r, 1) 
in terms of r (= rank of su{r, 1)) bosons. Whereas the comparison with the r + 1 boson 
realisation gives natural interpretation of various quantities and concepts. By explicit Lie 
algebraic calculation which goes quite parallel with probability theory, it is shown that the 
negative multinomial states are su{r, 1) coherent states belonging to discrete symmetric rep- 
resentations. Section 5 is for summary and comments. Appendix A serves to give general 
background of the paper, relating probability theory, coherent states and Lie algebra theory 
by taking elementary examples such as the ordinary coherent states and the binomial states. 
Appendix B also provides some basic elements like quantum mechanical generation of coher- 
ent states. A collection of two level atoms is discussed. It gives a good physical example of 
the binomial states and at the same time it provides simple interpretation of the H-P reali- 
sations as well as the relationship with the ordinary coherent states. Appendix C gives the 
higher rank generalisation of the results of the previous two Appendices. Here we advocate 
a seemingly ill-recognised fact that the multinomial states are coherent states of SM(r + 1) 
belonging to the symmetric representations. We stress, here as in the main text, the inter- 
play of probability theory. Lie algebra theory and quantum mechanics exemplified in various 
coherent states. Appendix D gives a short explanation of the negative binomial distribution 
as a distribution of "waiting time" . We adopt such notation as to reveal the essential features 
underlying this subject which sometimes results in deviating from the conventional notation. 

2 Negative Binomial State 

Let us start with the negative binomial distribution (For an elementary introduction of the 
negative binomial distribution from probability theory see Appendix D. For more details. 



3 



see for example, Chap. VI of [IIT 



B-{v;M)=i^^+^ y^-il-vT, n = 0,l,..., (2.1) 
in which < r/^ < 1 and M is a positive integer. This can be rewritten as 

(1 - vY''B-im M) = (^~^^^ i-v'r, n = 0, 1, . . . , (2.2) 

and it is easy to see that the right hand side corresponds to the power series expansion of 
(1 — rf y^'^^ the negative binomial expansion. Thus the normahsation 

oo 

Y,B-{r^-M) = l (2.3) 

n=0 

is obvious. From this it is also easy to see that the negative binomial distribution (and later 
the negative binomial states) can be defined for any positive number M. In this case we 
have to interpret 

fM + n-l\ T{M + n) 

\ n ) V{M)n\ ■ ^ ' ' 

Let us introduce the ^negative binomial state^ (NBS) by taking a 'square root' of the 
negative binomial distribution (|2 . 1| ) . To be more precise, we follow the analogy Poisson 
distribution <S=^ coherent state (for details see Appendix A): 

n! 

in which a > 0. Namely we define NBS 



n=o vn! 



oo 



\ve^';M)- = il-r^r'^^Y.\\ 



/M + n - r 



n 



ive^'^TWn), (2.6) 



in which {\\n) |n = 0, 1, . . .} are the number states of an oscillator: 



[6,6t] = l, 6||0)=0, H = ^^||0). (2.7) 



(The reason for using a slightly unconventional notation ||?t,) will become clear in the next 
section.) Then the number distribution in the NBS is the negative binomial distribution 

I inlve^'; M)"!^ = (1 - vY' ^ ~ = M). (2.8) 



n 

n2 



The condition < < 1 is necessary for the NBS to be normalisable. In the next section we 
will have a geometrical interpretation of the same condition as characterising the parameter 
space (the Poincare disk) of the sm(1, 1) coherent states. 



Next let us rewrite ( |2.6| ) {rjc = rje 

M ^ JM(M + 1) ■ ■ - (M + n- 1) 
\Vc; M)' = (1 - Ir/cH^ E -T^ -iVcrib^rWO). (2.9) 



„=o ^■ 



This can be reexpressed in the exponential form 

\Vc;M)- = (1 - |r/cn^exp[r/c7/C+] ||0), (2.10) 

in which 



/C+ = 6VM + = VM + N-lbl (2.11) 
Here use is made of the following identity |T2| 

(6^^(Ar))"||0) = (6^)"c/(0)c/(l)---^(n- 1)||0), with g{N) = ^^WTN, N = b^b. (2.12) 

Eq.(plOD and (|]TT]) reveal the su{l, 1) structure of NBS since /C+ and its hermitian conjugate 



IC- = VM + ATfe = bVM + N-1 (2.13) 



generate the sm(1, 1) algebra via H-P [jTO[ realisation of the discrete irreducible representation 
Q with the Bargman index M/2 : 

[/C+,/C_] = -2/Co, [/Co,/C±] = ±/C±, /Co = iV + ^, (2.14) 

and the "vacuum" ||0) is the lowest weight state: 

/C_||0) = 0, /CollO) = yllO). (2.15) 

It is easy to see tat ( |2.1(j| ) is expressed in the displacement operator form by using the 
disentangling theorem for sm(1, 1) : 

\r]e'^; M)" = exp [(c IC+ - QIC] ||0), Cc = e^^arctanhr/. (2.16) 

In other words the negative binomial states are su{l,l) coherent states in the definition 
of 0, 1^, 1^, although the sm(1, 1) structure is not obvious in the original definition of the 
binomial state ( p.6|) . It should be remarked that in contrast to the binomial states which 
cover all the coherent states of su{2) the negative binomial states give only part of the su{l, 1) 
coherent states. (There are other types of su{l,l) coherent states: for example those which 
are eigenstates of K_ 



^The generalised case with the real non-integer M gives the continuous irreducible representation of the 
universal covering group of SU{1, 1). 
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It should be remarked that the generating function of the negative binomial state 

oo 

G-iv;M;t) = 5:r5-(r/;M), |t| < 1 

n=0 



has a succinct "quantum" definition 



(2.17) 



G-{7]; M; t) = (r/e*'; M||t^||r/e'^; M). (2.18) 



As is well known in probability theory ||TT| the generating function is quite useful for calcu- 
lating various statistical quantities of the negative binomial states 0. 

3 Two Boson Formulation of NBS 

As with the binomial states discussed in Appendix A and B, the simplest way to understand 
the negative binomial states algebraically is to introduce two bosonic oscillators to express 
the sti(l, 1) generators as bilinear forms rather than the inhomogeneous forms as given ( |2.11| ), 
( |2.13| ) and ( |2.14| ). (We choose the formalism that the oscillators define the ordinary positive 



definite Hilbert space but the generators of the algebra reflect the non-compactness.) Let us 
introduce two bosonic oscillators 

[«i)4] = '^ifc, j,A; = 0,l, (3.1) 

and the Fock space 

fno fni 

a,|0,0) = 0, J = 0,1, |no,r^l) = ^^4^|0'0)• (3-2) 

Define 

K+ = ala\, K_ = aoai, Kq = -{Nq + + 1), A^^- = a]aj, (3.3) 
which satisfy su(l, 1) algebra 

[K+, K_] = -2Ko, [Ko, K^] = ±K^. (3.4) 

These operators either increase {K^) or decrease (K^) hq and rii simultaneously by 1 or keep 
them unchanged (Kq). In other words the above Fock space gives a reducible representation 
of su{l, 1) since the subspaces with different A = rig — are always separated. So we can 
restrict it as in the case of the binomial states 

A = no-ni = M-1 >0, (3.5) 
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in which M is a positive integer^^. Thus we arrive at the discrete representation of su{l, 1) 
with Bargman index M/2 

\M + n-l,n), n = 0,l,..., (3.6) 

with the lowest weight state 

M 

K^\M- 1,0) = 0, i^olM- 1,0) = y|M- 1,0). (3.7) 

Obviously this representation is irreducible. Since these states are uniquely specified by 
n = -n-i, we can identify them with the number states defined in the previous section ( |2.7| ) 
together with the "reduced" oscillator b and b^: 

\\n) = \M + n-l,n), n = 0, 1, . . . . (3.8) 

Thus we obtain the H-P representation of K± and Kq: 

K+^IC+ = bWM + N, K_ ^ /C_ = VM + iV6, K^^ 1Cq = N + (3.9) 

One advantage of the H-P type realisation as above is that it admits the generalisation to 
the continuous representation for non-integer M. 

The other group theoretical aspects of the negative binomial states are about the same 
as those in the binomial states. The physical and statistical properties of the NBS as well 
as their dynamical generation are discussed in some detail in our recent publication 0. The 
content of this section, though known in Lie algebra theory, can be considered to provide 
some Lie algebraic backgrounds for the probability distribution, which are new to the best 
of our knowledge. 

4 Negative Multinomial States 

The negative multinomial distribution is 

M-,{V ; M) = (1 - r7^)^^ ^^ + ^^^,|^•:+|;-"'^ ^ ■ ■ ■ r,^, (4.1) 
in which M is a positive integer and 

n = (no,ni, . . . ,nr), n' = (ni, . . . , n^) 77 = (?7i, . . . , r^^) G R'', 

< rf = ril^ Vril<l, n'! = rii! ■ ■ ■ ra^!. (4.2) 

^For the other sign of r.h.s. we can change the role of and 1. 

^In terms of the "waiting time" interpretation of the negative binomial distribution (see Appendix D) tiq 
is the total number of trials but the last (which is always a success), rii is the number of failures and M is 
the preset number of successes to be achieved. 
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In particular, the negative trinomial distribution reads 

M- ; M) = (1 - "IJ^;,;'!;)," '^;-'^'- («) 

This can be easily obtained from the negative binomial distribution 

B-{m M) = + " - - ^r', n = 0, 1, . . . , 



n 



by a binomial expansion 

2 2,2 2n 2ni 2n2 

and collecting appropriate terms. By repeating the same thing or by applying a multinomial 
expansion we arrive at the general form of the negative multinomial distribution ( ^.1| ). As 



we will see later this procedure also explains the generation of negative multinomial states. 

The negative multinomial state (NMS) is defined by taking a "square root" of the negative 
multinomial distribution (14.11): 



^'^ + "^ + --- + "^--^(.icr---(..cri|n'), (4.4) 



n'!(M-l)! 

in which the "reduced" states {||n') = \\ni, . . . , Ur) = 0, 1, . . .} are the number states of r 
bosonic oscillators: 

[b„bl] = 5,k, b,\\0)=0, j = l,...,r ||0) = |0,...,0), 

in') = ^^iio), {hr = br---br. (4.5) 

Vn'! 



It should be remarked that both negative multinomial distribution (|4.1j ) and state ( [4.4|) are 
also well defined for M positive real number. 

In order to show that the negative multinomial states are the coherent states of su{r, 1), 
we need to realise the algebra. Let us first construct su{r, 1) generators on the Fock space 
generated by r + 1 bosonic oscillators: 



[aj,al] = 6jk, aj|O) = 0, j = 0,l,...,r, |0) = |0, 0, . . . , 0), 

|n) = ^^=^|0), (a^)" = a[ ■ ■ ■ al . (4.6) 



Let us define the u{r,l) generators as bilinears in aj and a\: 

K+j = ala], K^k = aoOfc, l<j,k<r, 

Kjk = a]ak U^k^ 0), Nj = a]aj. (4.7) 
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It is easy to see that they leave the combination 

A = No-{Ni + --- + Nr) 

invariant and the above Fock space ( [4.6| ) is a disjoint sum of subspaces characterised by the 
value of A. As before, let us impose a constraint 

A = A/'o - (A^i + • • • + iV,) = M - 1 > 0, (4.8) 

in which M is a positive integer. Then for fixed M the restricted space provides an irreducible 
representation of u{r, 1). It has a lowest weight vector 

K_,|M-1,0,...,0) = 0, j = l,...,r (4.9) 

which is invariant under su{r): 

Kjk\M -1,0,...,0) = Nj\M -1,0,...,0) = 0, l<j,k<r. (4.10) 

Let us connect the states ||n') and |n). Each state in the above representation is uniquely 
specified by n' only and we identify 

||n') = |n) = |no, ni, . . . , n^), Uq = M + rii + ■ ■ ■ + rir — I 

or 

\\ni,...,nr) = \M + ni-\ h n,. - 1, ni, . . . , n^), and ||0) = |M - 1, 0, . . . , 0). (4.11) 

On these states the su{r, 1) generators are expressed inhomogeneously: 



J' 



K^. ^ JC^^ = byM + N^ + --- + Nr-l, /C+, = ^M + N^ + --- + Nr-lb 

Kjk^lCjk = h]hk, Afj = N,. (4.12) 

Note that the invariant subalgebra u{r) is expressed bilinear ly. 

It is not difficult to generate the negative multinomial states explicitly by applying the 
SU{r, 1) operator on the lowest weight state ||0). For simplicity and concreteness let us show 
this for the su{2, 1) case: 



|,,e-S,.e-^M)- = (l-^^)f f . ^^'t":rM"^7v'^ ^"^^0"(-^"^^^)"K--^)- 

nu^=o\ ni\n2l{M-iy. V / V / 

(4.13) 
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For the ^negative trinomial state' we only have to use the sm(1, 1) and su{2) disentanghng 
theorems for the two subalgebras spanned by (0, 1) and (1, 2) oscillators, respectively: 



= exp[e* ^ tanh r /Cio] exp [log(l - tanh^ r){Ni + M/2)] exp[— e * Hanhr/Coi], 

= exp[e^^2 tanr'/Csi] exp[-^ log(l + tan^ r')(A^i - A^2)] exp[-e-*^2 tanr'/Cis]. 
First let us choose r such that 

tanh^ r = rjf + rjl, 

to obtain 



oo 

M 



gKe-i/Cio-e--l/Coi)||0,0) = (1-7^2 _ ^2)f ^ 



(M + n- 1)! 



e^'V^f + %')"lKO), (4.16) 



„t^o^ (M-l)!n! 

which is a negative binomial state in the (0, 1) subspace. Next let us choose r' such that 



tanV = 



to obtain 



Thus we have shown that the negative trinomial states are the coherent states of su{2, 1) 
belonging to discrete symmetric representations. Note the parallelism with the negative 
trinomial distribution at the beginning of this section. The generalisation to higher rank 
cases is straightforward. One has to apply first su{l, 1) disentangling theorem and su{2) 
disentangling theorems in the following sequence of su{2) algebras spanned by (1,2), (2,3), 
. . ., (r — 1, r) oscillators. 

Before concluding this section let us remark that the generalisation of the discussion 
( |4.13| )-( p.l7| ) to the coherent states of su{r, s) is rather straightforward. 



5 Summary and Comments 

Stimulated by the well known fact that the photon number distribution of the ordinary 
coherent states is Poissonian, we have constructed quantum mechanical states which have 
the other well known probability distributions such as the binomial, multinomial, negative 
binomial and negative multinomial distributions as their particle number distributions. They 



10 



turn out to be the coherent states of the well known Lie algebras of su{2), su{r + 1), su{l, 1) 
and su{r, 1), respectively, belonging to certain symmetric representations. Interpretation of 
these coherent states in terms of probability theory is obtained and it is quite useful. At 
the same time Lie algebraic structure of these most fundamental probability distributions is 
revealed. 

The results of the present paper provoke many questions, to most of which we do not 
have answers yet. For example: What about the coherent states of su{r + 1) {su{r, 1)) 
belonging to the other representations? Are they also characterised by some probability 
distributions? The same question for the other Lie algebras, in particular the exceptional 
ones. Or the affine Lie algebras and other infinite dimensional algebras like the Virasoro and 
the w-algebras. . . ? 
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Appendix A Binomial States 

In the Appendix A-C, we reformulate the mathematical theory of coherent states for su{2) 
and su{r + 1) algebras. Most of the results are known in one way or another but we believe 
that the elementary exposition and the resulting explicit and concrete formulas and the 
emphasis on the connection with probability distributions are helpful and useful for most 
readers. It is also hoped that the comparison and the contrast with the compact cases will 
provide deeper understanding of the non-compact cases treated in the main sections. 
We follow the schematic path 




(A.l) 



by imitating the well known example of the Poisson distribution 



P„(a) = e-"' — , n = 0,1,2,.... 
nl 



(A.2) 



''Some coherent states can be constructed in terms of fermion oscillators 
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Let us introduce the "probability amplitude" by taking its 'square root'Q 

|«e*^) = e-'/2Vi^^|n), (A.3) 
n=o vn\ 

in which {|n)|n = 0, 1, . . .} are the number states of the ordinary oscillator 

[a,a^ = l, a\0) = 0, \n) = ^jL\0). (A.4) 



The origin of the additional phase factor e*^ is obvious, = ae'^^ae By using the last 
formula of ( |A.4| ) we can rewrite ( [A .31 ) as {ac = ae*^) 



\ac) = e^l"^l'/2 2 ^"^^"^^" [0) = e-l"^l'/2e°c^"'|0) = exp [aca^ - a^a] |0). (A.5) 

n=0 ^• 

At the last step use is made of the Baker- Campbell-Hausdorff formula. Eq.( |A.5| ) tells that 
the parameter space is the ordinary complex plane C, which is a coset space 

C = Heisenberg-Weyl Group/f/(l), (A. 6) 

in which the Heisenberg-Weyl Group is generated by a, and the identity operator. The 
stability subgroup U{1) is just the group of complex numbers of unit modulus, U{1) = 
{e'^ \9 : real}. 

The binomial distribution 

Br.{v;M)=(^^y'-il-vY-, n = 0,l,...,M, (A.7) 

is a well known elementary probability distribution related with binomial expansion^ 

1 = 1^ = (1 - a; + a;)^ = V (^]x^il - x)^"". (A.8) 

This gives the probability of n 'successes' among M times repeated Bernoulli's trials with 
the success probability < < 1. The associated 'probability amplitude' is 



M 



n=0 \ \ 



{7]e'yil~r]^)^\\n), (A.9) 



n 



in which {\\n) \n = 0,1, . . .} are the number states of the oscillator 



[b,b^] = l, b\\0) = 0, H = i^||0). (A.IO) 

Vnl 



""Apart from the well known example Klein-Gordon eq. Dirac eq., let us mention that the creation and 
annihilation operators , a are also 'square roots' of the oscillator Hamiltonian. In all these cases, including 
the probability amplitudes, the 'square roots' are complex in spite of the reality (hermiticity) of the original 
objects. 

Note that e" Pn{a) — is obtained by a power series expansion of e" . 
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(The reason for using a different oscillator b, from the above coherent state one a, a'^ and 
the slightly unconventional notation \\n) will become clear in Appendix B.) Let us call the 
state ( |A.9| ) the 'binomial state^ (BS) [1^. At first glance one might be tempted to give 

a phase to the second factor 

(V'l-r/V^^)^^-". (A.ll) 

But this is unnecessary since it is decomposed to an overall phase e*^^^ (which is immaterial) 
and e~*"^^ which can be absorbed by the redefinition of 6, 6 ^ 6 — 62. 
Next let us rewrite (|A.9|) {rjc = rje^^) 



\r]c; M) = (1 - \r]cr) 2 ^ 1 



Vc 



n=0 



(6t)"||0). (A.12) 



Then, by making use of the following identity [T^\ 

{b^giN)r\\0) = ib^rg{0)g{l) ■■■gin- 1)||0), with g{N) ^ v^M^, N = b% (A.13) 
we can write ([A.12|) in the exponential form 



\r]c;M) = (1 - \r]c\ ) 2 exp 



Vc 



1 - Ivcl 



l|0). 



Here 



J, = bWM - N = VM -N + lb^ 



together with its hermitian conjugate 

= VM-Nb = b VM - + 1 

generate the su{2) algebra via H-P [|10] realisation in the spin M/2 representation: 

M 

[J+, = 2 Jo, [Jo, J±] = ±J±, Jo = N-—, 



and the 'vacuum' ||0) is the lowest weight state 



M, 



J_||0) = 0, Jo||0) = -y||0). 



(A.14) 



(A.15) 



(A.16) 



(A. 17) 



(A.18) 



By using the disentangling theorem for su{2) we can rewrite (|A.14|) as 



lie''; M) = exp [Cc J+ - Cc J A l|0), Cc = e'^ arctan 



(A.19) 



In other words the binomial states are su{2) coherent states in the definition of 0, 
although the s-u(2) structure is not obvious in the original definition of the binomial state 
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( |A.9| ). Since all the irreducible representations of su{2) are exhausted by the representations 
( |A.1CI| )-( [A.18| ) for all non-negative integer values of M, the binomial states give all the su{2) 



coherent states. 

Before closing this Appendix, let us recall the fact that the binomial distribution tends 
to the Poisson distribution in a certain limit. Let M — > oo, — in such a way that the 
average value (n) is fixed: (n) = yfM = o? . Then for finite n 

OL 1 

5„(r/;M)^— e-° = P„(«). (A.20) 
n\ 

(Of course there are many other ways of showing this, e.g. in terms of the generating functions 
of these distributions.) We have the corresponding limit at the level of the "probability 
amphtude" (|A]|) 

\r]e''; M) ^ e-"'/^ £ (A.21) 

n=o yn\ 

namely the binomial state tends to the ordinary coherent state. This limit can also be 
visualised as a contraction of su{2) (|A.15| ),( |A.16| ) into the Heisenberg-Weyl algebra : 

riJ+ ah\ rjj^ ah. (A.22) 

Thus (pTTOD tends to 

\r]e'^] M) exp [ae'^ - ae"'^ b] ||0). (A.23) 

In Appendix B we will discuss the physical problem of dynamical generation of BS starting 
from certain Hamiltonian. This, in turn, will provide a mathematical framework in which 
(i) 5^(2) structure is more visible and, (ii) generalisation to the coherent states of su{r + 1) 
algebra, the ^multinomial states\ is straightforward. 

Appendix B Binomial States: Two Boson Formulation 

In order to discuss the generation of the binomial states, let us recapitulate the process of the 
physical generation of the ordinary coherent states, for comparison. This is an oversimplified 
model retaining only the most essential features of the coherent states. We focus on one 
particular mode of the photon since the system is decomposed into a sum of such subsystems: 

H = Ho + Hi, HQ = uja)a, Hi = j{t) + j{ty a, (B.l) 

in which a^,a are the creation and annihilation operators of the photon and j{t) is the 
classical current (with complex phase). The state vector in the interaction picture |'?/'(t))/ 
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obeys the equation of motion 

Hiit) = e'^^'H^e-'^^' = ]{t)e'^'a) +3{tYe-'^'a. (B.2) 
Let us suppose that the system is in the 'vacuum' |0) at t = 0. Then we obtain 

\i){t))i = Te-'/o«i(*')^*'|0), 

= gi!^{i)g(a(t)at-Q(i)*a)|Q^^ 



ait) = -% I jUy'^^^'dt', (B.3) 
Jo 

in which T is the time-ordering operator and Q{t) is a calculable function giving the imma- 
terial overall phase. 

For the binomial states let us consider a slightly different model consisting of a number 
of identical two level atoms (bosons). Let us also assume that the space extension of the 
system is not big compared with the wavelength of the photon corresponding to the energy 
gap and that the interactions between different atoms are negligible. The system can be 
described by the "spin" operators 

Hb = Hbo + Hbi, 

AI 

Hbo = ^^'^o\ 

M 

in which M is the number of the two level atoms. As is well known |16| a collection of identical 



particles can also be described by oscillators corresponding to each energy eigenstate. Let 
us denote the lower (upper) state and the corresponding oscillator by (1) [|: 



,j,al] = 6jk, j,k = 0,1. (B.5) 
Then the above Hamiltonian is equivalent to 



H'b = H'bo + H'b 

H'bi = ^,{t)a\aQ + fi(t)*alai, (B.6) 



— ^BO~^^Bl^ 

H'sQ = uiNi + ujqNo, Nj = a\aj, j = 0, 1, ui - uq 



^ In quantum optics situations one may call 'wiggler' photon and 1 'laser' photon p 
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and its Fock space is 

1^0, rii), no + ni = M. 
Now the su{2) structure is obvious, since 

J+ = alao, J- = ajai, Jq = ^(A^i - Nq) 



(B.7) 



(B.8) 



generate an su{2) algebra and the Fock space ( |B.7|) gives the M + 1 dimensional (spin M/2) 
irreducible representation corresponding to the Young diagram 



:□ M boxes 



with the lowest weight state 



M. 



J_|M,0) = 0, Jo|M,0) = -y|Af,0). 
The state vector in the interaction picture \ilj{t))j obeys the equation of motion 

n'Biit)\m)i, 

e'^Bo'H's^e-'^Bo' = fi{t)e''' J+ + /i(t)*e"^^* J_. 



(B.9) 



(B.IO) 



(B.ii: 



Let us suppose that the system is in the lowest weight state |M, 0) at t = 0. Then we obtain 



T exp 
Texp 



|M,0), 



|M,0). (B.12) 



Jo Jo 

Since g = T e""" Jo'^Bii^ y^ ^ SU{2) it can always be decomposed into g = exp(^J+ — 
C*J^) exp(zz/Jo) (C £ C, z/ G R) and the obtained state is the binomial state. For illustration 
purpose let us choose a special form of fi(t): 



flit) = ie 



et+ie 



(B.13) 



Then we obtain flj, |17 



\^{t))i = exp 
oc exp 



e*^tan(r7t) J J \M, 0) oc | tan(r7t)e'^; M). 



(B.14) 



However this is not exactly the same as the binomial state ( |A.14D given in the previous 
Appendix. In order to relate these two forms let us note that the state in the Fock space 
( |B.7| ) is uniquely specified by rii = -n, only: 



\\n) = \M — n,n). 



(B.15) 
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Let us understand that the states {\\n), n = 0, 1, . . . , M} are generated by the "reduced" 
single boson operators 6^ and b as in ( [A.10| ). Then the su{2) operators are expressed in terms 
of b'^ and b as 



•^+11^) = J+\M — n,n) = a[ao\M — n,n) , 



= Vn + WM -n\M — n — 1, 72 + — n ||n + 1), 

= b^yjM -N\\n), N = b^b, (B.16) 

namely 

J+ = bWM-N, = VM-Nb, Jo = N-^. (B.17) 

Thus we have naturally "derived" the H-P realisation of su{2) used in the previous Appendix. 
The lowest weight state in this notation is 

||0) = |M,0), X||0) = 0, JoiO) = -y ||0). 

At the end of the previous Appendix we have shown that the binomial state tends to 
the ordinary coherent state in a certain limit. Here we will show a result in an opposite 
direction. That is, the binomial states can be obtained from the ordinary coherent states 
with two degrees of freedom by appropriate 'slicing' or restriction. This reveals some features 
of the binomial states quite naturally. As before let us start with the corresponding result 
in the probability theory, which is rather elementary. A double Poisson distribution is given 
by 

fno,ni(ao,ai) = fn(Q:) = e~"°""i ° , \ , n=(no,ni), Q: = (ao,ai)- (B.18) 

no!ni! 

If we restrict it to a line 

no + rii = M, 

we obtain the binomial distribution up to normalisation: 



PM-n,n(ao, ai) = e -r^. 



(M-n)!n!' ' J^T^^' 



oc (^\''^{l-r]y'-'^ = B^{r];M). (B.19) 



n 



The same proposition at the level of the "probability amplitude" including the normali- 
sation can be easily obtained by considering the projection operator onto the representation 
space {^) [0: 

Pm= Y1 |no,ni)(no,ni|. (B.20) 
17 



With the aid of the resolution of unity (over-completeness relation) 



|"c)(Q!c| = 1, = («oc,«ic), d ajc = d{ajc)Rd{ajc)i, j = 0, 1 



(B.21) 



for the double coherent state 



|ac) = e' 2^ ^-p-r \no,ni) 

no,ni V^oini! 



we have 



M 



d'^aocd'^aic 



ac){occ\ Pm- 



By a change of variables 



c 



we have 



r.h.s. of (B.22) = J 



1 



(B.22) 
(B.23) 



00 



mo+mi 



C 



X E (^C'^l 



V^o!^! + led 
d\c (M+1)! 1 

vr (1 + iecP)*^+^ „,i;:=M VWW 



X ^ (ric^i 

no+ni=M 



(B.24) 



in which 



\ic) 



Y M 

1 A/ 



> n J 



. n , 



and the measure is 



(M + l)! d% 



c 



c) 



Ml 7r(l + |ecP) 



(B.25) 



(B.26) 



By introducing a parameter rjc = icl y 1 + l^cP? we can identify j.^c') as the binomial state 
|?7c; M) ( [A.9|) . This process shows elementarily that the parameter space of the binomial 
states is SU{2)/U{1) = CP^ {^c = cuc/aoc) obtained from {ctc=iceoc, aic)) by inte- 
grating out the overall factor (c- 
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Appendix C Multinomial States 

The multinomial distribution is 

M„(r7 ; M) = ^r^^ ■ ■ ■ //^"^(l - ry^)"", (C.l) 

in which 

n = (no,ni, . . . ,n^), no + H h = M, = (771, . . . , 77^) e R'', 

< rj^ = r]f + --- + r]'^ <1, n\ = UqW ■ ■ ■ rirl. (C.2) 
Let us first define the ^multinomial state' in the linear representation form 



\ric ; M) = T^^^iVicr ■ ■ ■ iVrcni - \ric\'r^», (C.3) 
in which the Fock states 

|n) = |no,ni, . . . ,?T,r), Uq + ni -\ h = M, (C.4) 

are generated by r + 1 bosonic oscillators 



75 "I] 



a,-, a 



S,k, %|0) = 0, j=0,l,...,r, |0) = |0,0,...,0), 



|n) = — |0), (a') = cii ■■■a|, . (<-^-5) 



n! 

Obviously the above Fock space ( |C.4| ) provides an irreducible representation of ^^(r + 1) 
with generators 

Jjk = a]ak, j 7^ k, Nj = a]aj, (C.6) 

in which Jj^ {j > k) are considered as shift-up operators. It is a symmetric representation 
corresponding to the same Young diagram as before: 

□□□•••m M boxes, (C.7) 

and the lowest weight state is 

|0') = |M,0,...,0), Jok\O')=0, J,fc|O')=0, j,k>0. (C.8) 

The last equation shows that the lowest weight state |0') is invariant under u{r). The 
dimension of the above irreducible representation is 

which is the same as the number of terms in the multinomial expansion, the number of 
the partitions of M into r + 1 non-negative integers and the number of M-th order partial 
derivatives of analytic functions of r + 1 variables. It should be remarked that there are 
other types of coherent states of su{r + 1) (r > 2) algebra belonging to the Young diagrams 
other than those given above ( p.7| ). They cannot be constructed by bosons only. 
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It is not difficult to generate tfie multinomial states explicitly by applying the SU{r + 1) 
operator on the lowest weight (energy) state 1 0') . For simplicity and concreteness let us show 
this for the su{S) case: 



M 



ni^=oV ^o!ni!n2! 

This process is essentially the same as the generation of negative trinomial state given in 
section 4. For the 'trinomial state^ we only have to use the su{2) disentangling theorems 
twice for two su{2) subalgebras spanned by (0, 1) and (1,2) oscillators: 

gr(e'«iJio-e-'»iJoi) -^-^^ 

= exp[e^^i tan r Jio] exp[log(l + tan^ r) (A^i + N2/2 - M/2)] exp[-e~^^i tan r Jqi] , 

gr'(e*^2J2i-e-'«2Ji2) -|^2) 

= exp[e*^2 i-j^nr' J21] exp[-^ log(l + tan^ r'){Ni - N2)] exp[-e-^^2 ^-^nr' J12]. 

We choose r and r' such that 

, 2 ^1 + , 2 / ^1 

tan r = ^ ^, tan r = 

1 - ^1 - ^2 ^1 
The generalisation to higher rank cases is straightforward. One has to apply sm(2) disentan- 
gling theorems in the following sequence of su{2) algebras spanned by (0,1), (1,2), (2,3), 
. . ., (r — 1, r) oscillators. 

To obtain the su{r + 1) multinomial states from the R + 1-fold coherent states is also 
straightforward. One only needs to develop clever notation to express the essential features 
succinctly. 

A multiple Poisson distribution is given by 



. _2no^2ni _ _ _ ^rir 

Pn(oc)^e~°'o 1 p— , n= (no,ni,...,n^), a = (ao, ai, . . . , a^). 

no! • • • riri 

(C.13) 



If we restrict it to a hyperplane 

no + rii + ■ ■ • + rir — M, or no — M — rii — ■ • ■ — rir, 
we obtain the multinomial distribution up to normalisation : 

r 

V c 



'Op- j=o 



P„(a) cx ^^f^---r?^'-(l-r/2)"o = M„(r7;M). (C.14) 
n! 
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The same proposition at the level of the "probability amplitude" including the normali- 
sation can be easily obtained by considering the projection operator onto the representation 
space (^) [0: 



Pm= Y1 \nQ,ni, . . . ,nr){no,ni, . . . ,nr 

riQ+niH \-nr=M 



(C.15) 



With the aid of the resolution of unity (over-completeness relation) 



r+1 



OLc){oLc\ = 1, etc = {aoc,aic, . . . ,arc) 
(fajc = d{ajc)Rd{ajc)i, j = 0, . . . , r 



(C.16) 



for the multiple coherent states 



lac) = e 



-|aclV2 



E 



(aoc)"" ■ ■ ■ {arcY 



no,...,nr. 



■|?2o,^l, ...,nr) 



we have 



M 



11^=0 d'^ajc 



71 



r+1 



M- 



By a change of variables 



faoc\ 
\arc ' 



cl 



/ 1 \ 

UrcJ 



(C.17) 



(C.18) 



we have 
Pm = 



X 



oo 

E 

mo,mi,...,mr= 



TT 



mo+miH hrrir 



Vm! \Jl + 1^ 



X E {no,nu 

no+ni+...+nr=AI 

d'^jc (M + r)! 

|2\M+r+l 



M 



1 + ie 



lie ^2C ' ' ' ^rC 



C\ 



E 



r?io+mi+...+mi-=M VIH' 
1 



X 



(1 + I^C 

{no, ni, . . . , nr\^Cw^2'c " " " C 



CaCc?---Ccl^o,mi,...,m.) 



rC 



(C.19) 



in which 



1 



(1 + |^^|2)A.//2^V n! 
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T.\l^^2b---Cc\W) 



(C.20) 



and the measure is 



d^i{^c,^c) - M! 7r-(l + |^^|2)r+i- ^^-^1) 



By introducing parameters rjjc = ^jc/yl + I^cPj can identify as the multinomial 
state |t7(^; M) ( p.3|) . This process shows elementarily that the parameter space of the multi- 
nomial states is SU{r + l)/f/(l) x SU{r) = CP*" {C,jc = c^jc/c^oc) obtained from C'^^ (ckc) 
by integrating out the overall factor (c- 

A few words about the multiple coherent states limit of the multinomial states. For the 
multinomial state ( |C.20| ) we let M ^ oo and C,jc ~^ while keeping the 'average' fixed, 



$,jcM = a -(J to obtain for fixed n' 



n' 

n' 



Like in the case of the binomial states one can express the states and the su{r + 1) 
generators in the "reduced" notation using only r boson oscillators. This gives rise to the 
generalisation of the Holstein-Primakoff realisation. But as remarked above it is applicable 
only to the symmetric representations. Because of the constraint 

nQ + nl^ Vrir = M, 

the state |n) is uniquely specified by 

n' = (ni,n2, ...,nr) 

only. So we identify 

||n') = Irii, 722, • • • , ^r) = |M — ^ rij, ni, . . . , Ur) = |n), (C.23) 
and introduce r independent boson oscillators 

[bj,bl] = 6,k, b,\\0) = 0, j = l,...,r. 



which create the "reduced" states 

(bt)"' 
n'! 

Then we have 



n') = i-i-||0), (bt)"' = 6r---&r, n'\ = no\---nr\. (C.24) 



Jjo = byM-Ni Nr, Jo, = \/M-iVi Nrbj, Jjk = b]bk. (C.25) 

Note that the "vacuum" ||0) is the lowest weight state and it is invariant under su{r) which 
is expressed linearly: 

^o,l|0)=0, Jjk\\0)=0. (C.26) 
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Before closing this Appendix, let us remark on the dynamical generation of the multi- 
nomial states. This is essentially the same as that of the binomial states. Let us consider 
a collection (total number M) of identical r + 1-level atoms (bosons). It is assumed that 
the interactions among different atoms are negligibly small compared with the interactions 
within the same atoms among different energy levels. As before the system is described in 
terms of r + 1 bosonic oscillators and the Hamiltonian at the zero-th order approximation 
is quadratic in the oscillators keeping the total number of atoms fixed. In other words the 
Hamiltonian is a hermitian linear combination of the u{r + 1) generators given in ( |C.6| ). If 
we assume that the system is in the lowest energy (weight) state |M, 0, . . . , 0) at t = 0, then 
at time t it is 

e-^-^*|M,0,...,0), 

which is a multinomial state since the time evolution operator e"*-'^* is an element of U{r + 1) 
and the U{1) part and the SU{r) is immaterial when they act on |M, 0, . . . , 0). 

Appendix D Waiting Time : Negative Binomial Distri- 
bution 

For those who are not familiar with probability theory, we give here a simple example in which 



the negative binomial distribution occurs. We follow Feller's textbook Let us consider 
a succession of Bernoulli's trials each of which has the probability of failure < rj"^ < 1. We 
ask a question: How long it will take for the M-th success to turn up? Here M is a positive 
integer. Since M-th success comes not earlier than M-th try, we denote by i?,^(?7;M) the 
probability that the M-th success occurs at the trial number M + n, n > 0. This occurs, if 
and only if, among the M + n — 1 trials there are exactly n failures and the M + n-th trial 
results in success: so that 

i?-(,;M)=(^^+;-')(l-r,yV". 

For the very unlucky the waiting time (n) can be infinite. This corresponds to the fact that 
the irreducible unitary representations of non-compact algebras are infinite dimensional. 
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